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A small variation of the Taylor Method and 
periodic solutions of the 3-body problem 
Oscar Perdomo 

Abstract 

In this paper we define a small variation of the Taylor method and a formula 
for the global error of this new numerical method that allows us to keep track of 
the round-off error and does not require previous knowledge of the exact solu¬ 
tion. As an application we provide a rigorous proof of the construction/existence 
of a periodic solution of the three body problem. Some images of this periodic 
motion can be seen at https://www.youtube.com/watch?v=fSmQyeKcj5k 



Figure 0.1: This is the orbit of one of the bodies in the periodic solution that we 
have mathematically shown to exist. The three masses are 100, 100 and 200 and 
the period of this solution is T = 96.241... We are assuming that the gravitational 
constant is 1. 
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Introduction 


1.1 The numerical Part 


When we think about using a numerical method to do a mathematical proof, the first 
thing that comes to our mind is that we need to consider the error of the numerical 
method. Very soon we realize that the standard formula for the error is no very useful 
due to the fact that it assumes that all the basic operations are being made with 
no error and this is computationally very expensive. As an example, if we consider 
the initial value problem y' = y — ^y^ with y(0) = ^ and we want to estimate y(^) 
using the Euler Method with h = we see that, even though we are considering 
the computational cheapest numerical method and we are only doing 30 iterations, 
a regular Computer Algebra System (Mathematica 10 in this case) will need 8696.99 
seconds to do these 30 iterations most likely because the final answer is a rational 
number of the form 2 with p and q integers both with 2607760525 digits, more than 
2.6 billions digits. On the other hand, if we allow the Computer Algebra System 
to have round-off error in every operation involved in each iteration, it becomes 
challenging to keep track of the error because easily, each iteration may have a few 
dozens of operations. In our example above, each iteration has 3 operations: one 
raising to the square, one product and a difference. We will exploit the fact that 
most of the Computer Algebra Systems (CAS) can compute, with a mathematical 
precision, the two integers i and j such that i < ri-\-^fr 2 < j where ri and r 2 rational 
number and r 2 > 0. The new numerical method that we are proposing in this paper 
allows us to work all the time with mathematical precision without paying the price 
of dealing with numbers that have huge expression using integers. Our numerical 
method will do the operations in each iteration with a mathematical precision and 
then at the end of the iteration, it will find with a mathematical precision a rational 
number that approximates the output of the iteration within a distance H. In order 
to be able to use the method, we find a formula for the error between the real value 
of the solution of the ODE and the approximation given by the numerical method in 
terms of the two values: h, the desired value for the step, and H, the desired value 
for the rounding in each iteration. This is done in Theorem (2.3). In general, when 


we try to estimate the error of a numerical method, a problem that we face is that 
we need to have some a-priori bounds of the solution of the differential equation 


that we do not know. An important aspect of Theorem (2.3) is that it does not need 
previous knowledge of the exact solution. 
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1.2 Periodic solutions of the three body problem 


Poincare showed that the three body problem has a chaotic behavior that makes 
it difficult to solve. For this reason, it is not surprising that the only explicit solu¬ 
tions were discovered more than 240 years ago by Euler in 1765 and by Lagrange 
in 1772. Recall that we have a periodic solution when the values of the positions 
and the velocities of the three bodies, after some time T, agree with the values of 
the positions and velocities at t = 0. Usually, when the values of a solution after 
T are within a small distance e from the initial condition, this solution is called 
numerically periodic. There is an enormous amount of numerically periodic solu¬ 
tions, for example, in 1975 [H], H. Henon showed a family of numerically periodic 
planar solutions of the three body with e = 10“^^. Despite the abundance of nu¬ 
merically periodic solutions, the task of showing that numerically periodic solutions 
are periodic is a difficult one. In 2000, Chenciner and Montgomery [T] showed that 
a numerically solution found by Moore in 1993, jl], was indeed periodic. This ex¬ 
ample represents (to my knowledge) the first example of a periodic solution that 
has a numerical image associated with it and does not have an explicit formula. 
It is important to mention that proofs showing the existence of periodic solutions 
have been found before, for example, Meyer and Schmidtl993 [3]. These solutions 
are usually near a bifurcation point of a family of solutions described with explicit 
formulas. After the paper by Chenciner and Montgomery, there has been more 
proofs showing that some numerical solutions are periodic. Among them, we have 
papers by Terracini and Ferrario, K.C Chen and Simo and Kapitza and Gronchi. 
An excellent account of the work done in this direction so far can be found in the 
site http://montgomery.math.ucsc.edu/Nbdy.html 


In this paper we will show that the numerical periodic solution given by the initial 
condition explained in Figure 1.1 is periodic. 


To describe the motion, let us assume that t represents seconds. Keep in mind that 
the units of time, mass and distance, have been adjusted so that the gravitational 
constant is 1. For this periodic solution the motion starts with two bodies, each 
one with a mass of 100, separated 20 units and the third body with a mass of 200 
right in the middle. If we reference the body in the middle as body 1, then, we have 
that the initial position of the body 1 is (0, 0, 0) and the initial position of the other 
two bodies are (10,0,0) and (—10,0,0). Body 1 will always stay in the z-axis while 
the other two bodies will move around the z-axis. During the first to = 2.673... 
seconds. Body 1 moves straight up to the position (0, 0, 2.453 ...) = (0,0, zi), this is 
the farthest up this body will get. Simultaneously, the other two bodies move down 
and around the z-axis; after these first to seconds they are both zi units below the 
x-y plane, they are ri = 9.43 ... units away from the the z axis; and they have made 
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Figure 1.1: There is a periodic solution with period T = 96.241..., a = 4.3170 ... 
and b = 1.4903 .... 


a ^-rotation (70 degrees) with respect to the z-axis. At every instance, the distance 
to the z-axis of these two bodies is the same and it is twice the distance between 
them. In this way, after the first to seconds Body 2 is at (n cos ri sin —zi) 
and Body 3 is at (—ri cos —ri sin —zi). During the next to seconds the body 1 
moves down getting back to the origin and the other two bodies will simultaneously 
rotate, move up and move apart from the z axis; they will rotate another and 
they will go back to be 10 units apart from the z-axis. Notice that after 2to seconds 
the relative positions with respect to each other are the same as in the starting 
position. The only difference is that Body 1 is now going down and the other two 
bodies are going up. With respect to a fixed reference frame, after 2to seconds, the 
positions of all three bodies differ from the starting position by a rotation of 
The next to seconds Body 1 will reach its lowest point (0,0,—zi) while the other 
two bodies will reach the highest point rotating another Finally after another to 
seconds the relative positions with respect to each other are the same as the starting 
position, and this time the body 1 is going up as it was at the starting position. 
With respect to a fixed frame, the positions after 4to seconds differ from the starting 
position by a rotation of Doing 9 more of these 4to cycles will bring the three 
bodies to the starting position. In this process, bodies 2 and 3 will have completed 
7 rotations around the z-axis. The trajectory of Body 2 is shown in Figure 0.1 
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Each color represents a cycle of 4to seconds, for this reason there are 9 colors in the 
picture. Body 3 does not share the same trajectory as Body 2. The trajectory of 
Body 3 is the reflection with respect to the x-y plane of the trajectory of the body 
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We would like to emphasize that the goal of this part of the paper is not to show 
the existence of this type of solutions: This was done by Meyer and Schmidt [3], nor 
to show numerically periodic solutions of this type: We can see similar images of 
solutions like the one we are showing (when the three masses are the same), in the 
work by Yan and Ouyang [7j . The goal of this part of the paper is to give a rigorous 
proof that a particular numerically periodic solution is indeed periodic. 

We will reduce the proof of the periodicity of this solution to show that three func¬ 
tions defined in an open set of B? must vanish simultaneously. The variables in 
the domain of these three functions are given by triples (t, a, h) where a and b are 
explained in Figure o 

The author has found some numerical solutions of the same type, not only for the 
3 body problem but for the n-body problem. Some aspects of these solutions have 
been posted online: 

The link https://www.youtube.com/watch?v=PtEMb6Rvflg shows a periodic solu¬ 
tion of the 6 body problem. 

The link https://www.youtube.com/watch?v=2Wpv6vpOxXk shows a periodic so¬ 
lution of the four body problem. 

The link https://www.youtube.com/watch?v=hjQplP09560 shows a periodic solu¬ 
tion of the three body problem. 

The images in the videos were generated by solving the differential equation that 
governs the n body problem and then posting 10 pictures per second, in this way, 
the time in the video is proportional to the time t of the solution. 


The technique used in this paper to prove periodicity is different to the one used 
by Chenciner and Montgomery, they used variational methods. There are five in¬ 
gredients in our proof that the solution that we are considering is periodic: (i) 
The Poincare-Miranda theorem. Theorem 4.4, which is essentially a generalization 
of the intermediate value theorem, (ii) A symmetry result that allow to integrate 
the ordinary differential equation (ODE) over a quarter of a period instead of the 
whole period. Lemma 6.2 (iii) The Round Taylor Method to solve differential 


equations, that allow us to estimate the values of the functions that we are con¬ 
sidering. Section (iv) A lemma related with the the implicit function theorem 
that allow us to find a set where the solution of a system of equations of the form 


5 




fi{t,a,b) = 0, f 2 {t,a,b) = 0 is given by only points on a connected smooth curve, 
Section (v) A theorem that tell us how to compute the partial derivative of an 
ODE with respect to the initial conditions and parameters in the ODE. [2]. 

The author would like to express his gratitude to Andres Rivera, Richard Mont¬ 
gomery, Carles Simo and Marian Anton for their valuable comments on this work. 


2 The Round Taylor Method 

2.1 The Round Taylor Method 

We will be using a small modification of the Taylor method of order m to estimate 
solutions of differential equations. This modification will allow the computer that 
implements the method to work all the time with rational number and its square 
roots instead of working with approximations of numbers. In this way, the method 
is free of error coming from decimal approximation, this is, this method is round -off 
error free. Let us call this method the Round Taylor Method. Assume we have the 
ordinary differential equation 

Y'{t) = f{Y{t)), Y{0) = yo (2.1) 

where Y{t) = (Yfft),... ,Yn{t))'^ and f{y) = {fi{y), ■ ■ ■, fniy)V ■ Here v'^ stands 
for the transpose of the vector v. Given two positive rational number h and R, we 
define the sequence of points in that starts with zo = yo and follows using the 
recursive formulas 


2/1 

= zo +f{zo)h +Fffzo)—+ . 

w 

' • ^ m- 

hm 

- 1 (^ 0 )^ 

m\ 

Zl 

= a rational number such that 

\zi-yi\< H 

2/2 

= zi +f{zi)h +Fffzi)^ + . 

W 

• • m- 

hm 

m\ 

Zi 

= a rational number such that 

\zi -yi\<H 

Vi+i 

= Zif {zi)h + Fffzi)-^ + .. 

• Fm- 

hm 

m\ 


where Fi = Dff, F 2 = DFif,... Fm-i = DFm- 2 f- Here D denotes the derivative 
operator that takes a function from R" to R” to the matrix which columns are the 
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partial derivatives with respect to the variables. In order to define Zi from yi, we 
may use -and will use in this paper- the Floor function [rcj that assigns to a real 
number x the largest integer not greater than x. When x = (xi, ..., x^) is a vector, 
[xj = ([xij,..., [xfcj). It is clear that a posible choice for Zi is Zi = . 

As an example, when we use the Round Taylor method of order 1 for the differential 
equation given in the introduction, this is, n = 1, f{y) = y — yo = ^. When 
h = 1/100 and H = 10“®, we obtain that yo = 1/2 = zq and the next 10 values of 
{yi, Zi} are given by 


121 

252083 

240 

500000 

38127028661111 

12709 

75000000000000 

25000 

96109156319 

256291 

187500000000 

500000 

38762401423319 

16151 

75000000000000 

31250 

152668926449 

521109 

292968750000 

1000000 

52541490803373 

262707 

100000000000000 

500000 

13243698510717 

529747 

25000000000000 

1000000 

160232709115991 

534109 

300000000000000 

1000000 

161549754576119 

538499 

300000000000000 

1000000 

162875215826999 

542917 

300000000000000 

1000000 


Let us start the process of finding the formula for the error of this numerical method. 


Lemma 2.1. The sequence given by the recursive formula qi+i = {1 + L)qi+p (with 
L / 0/ satisfies that qt = (1 + L)^ — £. 


Proof. The proof follows by induction. Clearly the formula works for A: = 0. Now, 
assuming the formula works for k >0, we have: 


Qk+i — +LjQk+P — +L) { (1 +-^) fi (1 + -^) 
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This finishes the proof of the lemma. 


□ 


There are different norms that we can use in the set of matrices, in order to establish 
the the one that we are using, we state the following lemma. 

Lemma 2.2. Let A be an n x n matrix, if \A\‘^ denote the sum of the square of its 
entries, then for any veetor v in i?” we have that \Av\ < |vl||u|. 

Proof. If A^ denote the columns of A and Vi denote the entries of v, then we have 
that 


i=\ i=l 

□ 

Theorem 2.3. Let h and H be two positive numbers and k a positive integer. Let 
us consider the sequences {( 2 /i, by the Round Taylor Method associated 

with the values h and H for the ODE Y\t) = fiY (t)), Y (0) = yo, deseribed in the 
beginning of this seetion. Let us assume that we ean find constants Mq, Mi, ... Mn, 
Kq, ... ,Km-i and sets Ui = {{ui,...,Un) : bi < ui < ci,...,6„ < Un < Cn] 

and U 2 = {(iti, • • • ,Un) : bi — € < ui < ci + e,... ,bn — c < Un < Cn + e} with 

e > Moh + H, where H = - l) and M = M'l + • • • + M 2 and 

L = Ko + Ki^ + --- + If 

• The map f and all its partial derivative with order less than m + 2 are contin¬ 
uous in an open set that eontain the elosure 0 / 1 / 2 . 

• Zj G Ui for j = 0,...,k. 

• \Fmiiu)\ < Mi for i = 1,... ,n, for all u G U 2 . Here Fm = {Fmi, ■ ■ • ,Tmn)- 

• \Df{u)\ < Kq and \DFi{u)\ < Ki for i = 1,... ,m — 1 and u ^ U 2 . 

• \fi{u)\ < Mq for i = 1,... ,n and u £ U 2 . 

then, the solution Y(t) of the system of ordinary differential equation is defined on 
[0, hk] and for any positive integer j < k, we have that \Y(jh) — Zj\ < H. 



Proof. Let us start by checking that for any pair of points ui and U 2 in U 2 we 
have that \f{u 2 ) — < Ko\u 2 — ui\ and \Fi{u 2 ) — Fi{ui)\ < Ki\u 2 — ni| for 

i = 1,... m — 1. We have that 


\f{u2) - 


df{ui + t{u 2 - Ul)) 


dt 


dt 


/ Df {ui + t{u2 - Ui)) {U2 - Ui) dt 

Jo 


< Kq\u 2 - Ui 


The proof of the inequalities \Fi{u 2 ) — Fi{ui)\ < Ki\u 2 — wi|, for i = 1,..., m — 1 
is similar. For any non negative integer I < k such that Ih is in the domain of Y 
let us denote by xi = Y{lh). Notice that anytime Y{t) is in the closure of U 2 for 
all t G then it is possible to extend Y{f) to an interval of the form [0,ti + 5] 

with (5 > 0. We will show that for any I < k, Y{t) is defined on [0,lh], Y{t) G U 2 
for all t G [0,hl] and \Y{lh) — zi\ < H. Clearly the result hold for ^ = 0 because 
\xo — zq\ = 0. Let us assume that Y{t) is defined for all t G [0,jh] with Y{t) G U 2 
and \xi — zi\ < H for all 0 < I < j < k, we will prove the theorem by showing that 
Y(t) is defined for all t G [0, {j + l)h] with Y{t) G U 2 and |xj+i — Zj+i\ < H. Let us 
show that for any t G [j/i, (j + l)h], Y{t) G 112- Since Y{t) is continuous and U 2 is 
open, then Y{hj + r) G 1/2 for small positive values of r. Let us show by contraction 
that Y{hj + r) ^ U 2 for all r G If Y{hj + t) ^ U 2 for all r G [0,/i], the 

we can find r* such that Y{jh + t) ^ U 2 for all r G (0,r*) and Y{jh + r*) ^ U 2 . 
Therefore, writing Y = (Yi,..., Yn), we can find a positive integer w < n such that 
either Y^ijh + r*) = — e or Yw{jh + r*) = + e and Ye{jh + r) G {be — e, Cg + e) 

for all positive integers e < n and all r G (0, r*). Denoting Zj = ..., Zj^), we 

have that for some ^ G (0,r*), 


Yy,{jh + r*) - Zj^ 


< {YwiJh + T*) - Y^{jh)\ + \Y^{jh) - Zj^ 

< \YUhj + 0\r* + H 

< \f^{Y{jh + 0)\r* + H 

]\P(}h T FI € 


Since Zj G Ui, the inequality above contradicts the fact that either Yw{jh + r*) = 
bw — € or Yw{jh + r*) = + e. This contradiction shows that Y{t) G U 2 for all 

t G [0, (j + l)h], in particular Xj+i G U 2 . Let us prove now that \xj+i — 2 :^+ 1 1 < H. 
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For any i = 0,..., j' + 1 we have that 


l^j+i ■2'j+il ^ 

< 


< 


l^i+i 2/j+i| “1“ l2/*+i -^j+il 

h? /i™ 

\xi + f{xi)h + Fi{xi)— H - h Fm-i{xi) — j- 

2\ ml 


Xi — Zi\ {1 hL) + - -ywT7 + H 

[m + Ij! 


yi+i \ + 


(m + 1)! 


+ H 


Let us define p = 
can show that \xi 


+ H, qo = 0 and Qi+i = (1 + hL)qi + p. By induction we 
— Zi\ < qi for alH < j + 1. Using Lemma (2.1) we obtain that 


\Y{{j + l)h) - Zj+i\ = \xj+i - Zj+i\ < ^{{l + hiy^^-1) 

< :^(e“^-l) = # 

tllj 

This finishes the proof. 

□ 

The following theorem is well known. A reference for the particular case when / 
and g are Lipschitz real value functions can be found at Earl A. Coddington, An 
Introduction to ordinary differential equations -Dover Publication - 1989. A refer¬ 
ence for a more general case can be found at Herbert Amann, Ordinary Differential 
Equations An Introduction to Nonlinear Analysis - Walter de Gruyter - 1990. 

Theorem 2.4. Let us assume that f,g:U C —)• are functions defined on 
an open convex set U such that \Df{u)\ < K, \Dg{u)\ < K and \ f{u) — g{u)\ < e. 
If y{t) and z{t) satisfy y{t) = f{y{t)) and z{t) = g{z{t)) respectively, then 

\y{t) - zit)\ < | 2 /(to) - ^ - l) 


3 On the implicit function theorem 

Let us consider the set S = {x = (xi,X 2 ,X 3 ) : /i(x) = 0, f 2 {x) = 0} where 

/i and /2 are real value functions defined on an open set of R^. By the Implicit 
Function Theorem we know that if po £ ^ and V/i(po) x V/ 2 (po) is not the zero 
vector then there exists an open set U that contains po such that S n 1/ is given a 
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regular connected curve. The following lemma give us an estimate on how big this 
open set U can be, under the assumption that we know that there is point in S in 
a small box with dimensions 2 //i, 2 fj ,2 and 2 // 3 . 

Theorem 3.1. Let us assume that fi, f2'R^^R are smooth functions, ni, ^2, hs, 
ei; ^2, ^1, £2; <^27 <^1, 62 urc positive numbers such that for i = 1,2, < e* < 6 i, 

P^ = "-^-h3>0 where mi = and m2 = and 


U = {{x,y,z) : |x| < ei, \y\ < 62 , jz] < p] wherep < min{pi,p 2 ] ■ 

If there exists po = {xo,yo,zo) such that fi{po) = f 2 {po) = 0 with |xo| < pi, 
1 2 / 0 1 < P 2 , \zo\ < P 3 and for every p € U we have that 


5i < 


dxi 


< ^ 1 , 


dx2 


< ei, 


dx3 


< ei, 62 < 


df2 

dx2 


< h, 


< 9/2 

5x1 


< €2, 


^2 

dx 


then, for every zi with \zi\ < p, there exists a unique {xi,yi) such that (xi,yi, 2 ;i) G 
U is a solution of the equations fi{x,y,z) = 0, f 2 {x,y,z) = 0. 


Proof. Let us denote by .^ = (^ 1 , C 2 , ^s) = V/i x V/ 2 . Direct computations show that 
161 < €i (52 + e 2 ), 161 < e 2(6 + ei), and | 6 | > 66 -eie 2 - Let a{t) = {x{t),y{t), z{t)) 
be the integral curve of the vector field ^ such that a(0) = po- Since | 6 | > 0 we have 
that the vector field ^ never vanishes on U, therefore there exist r 2 > 0 and Ti < 0 
such that a{t) G U for all t G (Ti,T 2 ) and a(ri) and a{T 2 ) are in the boundary 
of U, in particular we have that either x(T 2 ) = ± 6 , 2 /( 6 ) = ±6 or z{T 2 ) = ±p. 
We will prove that 2 ;(T 2 ) must be either p or —p. Using the fact that x = fi{a{t)), 
y = 6 (ck(^)) and i = 6 (<a( 6 ) wo obtain that. 


X.. €1(62+ €2) , i2/| . 02(6 + ei) 

V < -= mi and v < - 

z 0102 - 6162 z O1O2 - £102 


m 2 


(3.1) 


Recall that z{t) is one to one because its derivative never vanishes, 
by h the inverse function of z. We have that h{z{t)) = t. If 51 (r) = 

£l(<l2+e2) 

-eie2 ■ 


gi{zii) = xo and \g'i{T)\ < Therefore, 


Let us denote 
x{h{T)), then 


x{T2) - xol = \g{z{T2)) - 51 (^ 0 )! = \9'ix*)\\z(T2) - zo\ < mi6(r2)| + mip3 
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and then |x(T 2 )| < //i + mipi +mi/i 3 = ei. Likewise we can show that \y(T 2 )\ < € 2 - 
Therefore we must have | 2 ;(T 2 )| = p. The same arguments show that | 2 :(Ti)| must 
be either p or —p. Since the function z{t) / 0 and [Tl,r 2 ] is connected then we get 
that for every zi with l^ij < p, there exists a pi = {xi,yi,zi) = a{ti) G U such that 
= 0 and / 2 (pi) = 0. Let us prove that pi is unique. If p 2 = {x 2 -,y 2 -,zi) G U 
and p 2 / pi, then either |x 2 — xi| < \y 2 — yi| / 0 or \y 2 — yi\ < |x 2 — xi| / 0. If 
|x 2 — xi| < \y 2 — yi\ / 0. Since we are assuming that / 2 (xi,X 2 ,X 3 ) = 0 then, 


l/2(x2,y2,^;i)| > |/2(x2,y2,^;i) - /2(x2,yi,^i)| - |/2(x2,yi,^;i) - f2{xi,yi,zi)\ 

> I 2/2 - yi |52 - |X 2 - Xi|e 2 > 0 

Likewise we can show that if \y 2 — yi| < |x 2 — xi| / 0, then |/i(p 2 )| > 0. Therefore 
the only solution of the system fi{p) = f 2 {p) = 0 on the set U are those points in 
the curve a : [Ti,T 2 ] —)• U. This finishes the proof. □ 

For any u = (ui, U 2 , U 3 ) and v = (ui, U 2 , U 3 ), let us denote by u-v = U 1 U 1 +U 2 X 2 + 143 X 3 
the dot product of u and v. Given a function f : ^ R, we denote D^fi = 'V fi-v. 

By doing an orthogonal change of coordinates we have the following corollary, 


4 The differential equation, a symmetry result and the 
Poincare Miranda Theorem 

Let us start introducing the differential equation for the subfamily of solutions of 
the three body that we are considering in this paper. 

Proposition 4 . 1 . If f, r and 9 satisfy the initial value system of differential equa¬ 
tions 


f 


400 ^ lOOa^ 


25 



(4.1) 


with r(0) = 10, /(O) = 0, r(0) = 0, /(O) = b, 6 { 0 ) = 0, where s = + 4/^ 

then, 
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x{t) = (0,0,/(i)) 

y{t) = (r(t)cos(6»(t)),r(t)sin(6»(t)),-/(t)) 

z{t) = {-r{t)cos{e{t)),-r{t)sm{9{t)),-f{t)) 


is a solution of the 3-Body problem with the mass for the body moving according to 
x{t) equal 200 and the masses of the other two bodies equal 100. We are assuming 
that the gravitational constant is 1. 

Proof. Notice that \x — y\ = |x — 2 ;| = h and \y — x\ = 2r, we can check that if / and 
r satisfy the ODE, then x{t) = 100 + 100|^f^, y{t) = 200|^f^ + 100|^f^ 

and z{t) = 200 + 100 |^~^^3 . For more details see |6] or [3] for example. □ 

Lemma 4.2. Let us assume that f and g satisfies the ordinary differential equation 

f = </>(/, 5 ) 
g = ^{f,9) 

with /(O) = 0, /(O) = a, g{0) = r > 0 and g{0) = 0 and f and ^ smooth functions. 
If(p{-f^9) = -4>{f,9) and C{-f,g) = f.{f,g) then, f{t) is odd and g{t) is even. 

Proof. Let us consider the functions f{t) = —f{—t) and g{t) = g{—t). A direct 
computation shows that /(O) = 0, /(O) = a, ^(0) = r and g(0) = 0. Moreover we 
have that 

Rt) = -fi-t) = -4>{f{-t),g{-t)) = -(j){-f{t),g{t)) = 4>{f{t),g{t)) 


and 


g{t) = g{-t) = f.{fi-t),g{-t)) = f.i-f{t),g{t)) = ^{f{t),g{t)) 

By the uniqueness of the solutions of ordinary differential equations we get that 
f{t) = f{t) = —f{—t) and g{t) = g{t) = g{—t). This finishes the proof. 

□ 
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Lemma 4.3. Let us assume that f and g satisfies the ordinary differential equation 

f = Hf^g) 
g = C{f,g) 


ff /(O) = 0 and ^(0) = 0 and (j) and ^ smooth functions, then, f{t) and g{t) are 
even. 

Proof. Let us consider the functions f{t) = f{—t) and g{t) = g{—t). A direct 
computation shows that /(O) = /(O), /(O) = 0, g(0) = (/(O) and ^(0) = 0. Moreover 
we have that 


Rt) = f\-t) = (p{f{-t),g{-t)) = f){f{t),g{t)) 

and 


g{t) = g{-t) = i{f{-t),g{-t)) = fif{t),g{t)) 

By the uniqueness of the solutions of ordinary differential equations we get that 
f{t) = f{t) = f{—t) and g{t) = g{t) = g{—t). This finishes the proof. 

□ 


This is the Poincare-Miranda theorem for two variables 

Theorem 4.4. Let U be an open set in B? that contains the rectangle [ai, 02 ] x [ti,t 2 ] ■ 
Let us further assume that F : U x I —)■ R and G : U x I —)■ R are continuous 
functions. If 

(i) F{a, ti) > 0 and F{a, t 2 ) < 0 for all a G [oi, 02 ], 

(a) G{ai,t) < 0 and G{a 2 ,t) > 0 for all t G [^ 1 ,^ 2 ]? 

then, there exists a point (ao,to) £ [ 01 , 02 ] x [^ 1 ,^ 2 ] such that F{ao,to) = 0 = 
G(ao,to)- 


5 The solution of the ODE as a function of the time 
and the parameters a and b. 

The main tool used to show the periodicity of the solution of the three body problem 
is to understand the functions in the solution of an ODE as functions of t, a and b. 
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We this in mind we define, 


Definition 1. We will denote by F{t, a, b) = f{t) and R{t, a, b) = r{t) and 0(t, a, b) = 
6{t) the solution of the system (f-D with initial conditions 


/(O) = 0, /(O) = b, r(0) = 10, r(0) = 0 and 0(0) = 0 

We will denote by R = ^ and F = ^ and in general if S is a function of the 
variable {t, a, b) then S = ^ 

It is well know that the function F and R has continuous partial derivatives [2] and 
they obey a differential equation. 

Theorem 5.1. If we denote by Fa = Ft, = ^, Ra = ^ and Rt, = ^, then 
for any fixed values a and b, 


Y{t) 


{F{t, a, b), F{t, a, b), R{t, a, b),R{t, a, b), Fa{t, a, b), Fa{t, a, b), 
Rait, a, b),Ra{t, a, b), Fb{t, a, b),Fb{t, a, b), Rb{t, a, b),Rb{t, a, b)) 


satisfies the differential equationY = P{Y) with P{xi,... ,xi 2 ) = P{x) = {x), ■ ■ ■, P^'^ (x)) 

where 
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P^{x) 

P‘^{x) 

P^{x) 

P^{x) 

P^{x) 

P^{x) 

P\x) 

P^{x) 

P^{x) 

pi°(x) 

P^^{x) 

P^^{x) 


X2 

400xi 

(4xf + x|) 3/2 

X4 

lOOa^ _ 25 _ 200x3 

x| x| (4xf + x|) 3/2 

X6 

3200x5xf 1200x3x7x1 400 x|x5 

(4xf + X3) 3/2 (4x2 + x|) 3/2 (4^2 -|- x|) 3/2 

X8 

300a2x7 200a 2400x1X3X5 50x7 600x|x7 200x7 

X3 X3 (4xf + x|) 3/2 X3 (4xf + x|) 3/2 (4^2 + x|) 3/2 

XlO 

3200x9X3 1200x3X11X1 400 x|x 9 

(4xf + x|) 3/2 (4^2 + x|) 3/2 (4^2 -|- x|) 3/2 

a;i2 

300a2xii 2400x1X3X9 50xii 600x|xii 200xii 

X3 (4xf + x|) 3/2 X3 (4xf + x|) 3/2 (4x| + x|) 3/2 


with initial conditions 


Y{0) 


( 0 , 6 , 10 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 0 , 0 ) 


6 Main result 

In this section we explain the arguments used in the proof of the periodicity of the 
solution and we will also show how to reduced this proof to the proof of 4 lemmas. 
With the intension of not cutting the flow of main ideas, we will prove these lemmas 
in a different section. 
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6.1 Special values used in the proof 


before we continue explaining the main result we would like to defined some con¬ 
stants. 

43170475352787 ^ _ 17 _ 1197 

“ 10000000000000 “ “ 50000000 “ 100000000 

1490359743 1 

° “ 1000000000 “ 50000 

13366894627923 ^ _ 1 _11 

“ 5000000000000 ^“ 2500000 ~ 2000000 

6.2 Main Theorem: Periodic solution 

Our main theorem shows that for some values of a, 6 and T of the form a = 
4.3170 ..., b = 1.4903 ... and T = 96.24 ..., the solution of the three body problem 
is periodic. For this periodic solution the orbit of the body with mass 200 goes up 
and down on the ^;-axis, and the other two bodies moves around the z-axis; the orbit 
of one of them is shown in Figure [Odj 

Theorem 6.1. There exist a triple {a,b,F) with |a —ao| < sa + da, \b — bQ\ < sb and 
\t — to| < Q{ds + dt), such that the solution of the three body problem given by 


x{t) = {0,0,F{t,a,b)) 

y{t) = {R{t,d,b) cos{Q{t,a,b)), R{t,d,b) sm{Q{t,d,b)),—F{t,d,b)) 

z{t) = {—R{t, a, b) cos(0(t, d, b)), —R{t, d, b) sin(0(t, d, b)), —F{t, a, b )), 

is period with period 36t. Moreover, we have that Q{d, b,t) = 

6.3 Symmetry lemmas and reduced periodic solutions 

A solution of the three body problem is called reduced periodic if the functions that 
provide the distances between the bodies are periodic with the same period Tq. In the 
subfamily of solutions that we are describing in this paper, reduced periodic means 
that for some Tq > 0, R{t + kTQ,a,b) = R{t,a,b) and f{t + kTo,a,b) = f{t,a,b) 
for any integer k. The following theorem makes easier the task of finding reduced 
periodic solutions. 
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Lemma 6.2. If for some ti > 0 we have that R{ti,a,b) = 0 = F{ti,a,b), then the 
functions f{t) = F{t, a, b) and r{t) = R{t, a, b) are periodic with period 4ti. 


Proof. By Lemma (4.3) we have that F{ti — t,a,b) = F{ti + t,a,b) and R{ti — 
t, a, b) = R{ti + ta, b). It follows that —F{ti — t, a, b) = F{ti + t, a, b) and —R{ti — 
t,a,b) = R{ti + t,a,b), therefore R{2ti,a,b) = R{0,a,b) = 0 and F{2ti,a,b) = 
F{0,a,b) = 0. By Lemma 4.2 we have that f{—t) = —f{t) and r{—t) = r{t). 
Using Lemma 4.2 one more time we obtain that f{t + 2ti) = —f{2ti — t) and 
r{2ti + t) = r{2ti — t). It follows that /(4ti) = /(O) = 0, /(4ti) = /(O) = b, 
r(4ti) = r(0) = 10 and r(4:ti) = —f(0) = 0. Since the values of the solutions / and 
r at t = 4ti agrees with those at t = 0, then the lemma follows. □ 


By the previous theorem, a point (ti, oi, bi) satisfying R(ti,ai, bi) = 0 = F{ti, ai,bi) 
defines a reduced periodic solution. If additionally 0(ti,ai,6i) = ^ with p and 
q > 0 integers, then the solution is not only reduced periodic but periodic with pe¬ 
riod 4gti. Moreover, by the implicit function theorem if the cross product between 
the gradients of the functions R and F at (ti,ai,5i) does not vanish, then there 
is a curve of points that solve the equation R{t, a,b) = 0 = F{t, a, b) and then, we 
obtain a family of reduced periodic solutions. In the case that the function 0 is not 
constant along this curve of points in the space that represent reduced periodic so¬ 
lutions, then we obtain infinitely many periodic solution due to the fact that on any 
open interval there are infinitely many numbers of the form with p and q integers. 
The Implicit function theorem also tell us that there exists an small open set around 
(ti,ai,6i) such that all the solutions of the equations R{t,a,b) = 0 = F{t,a,b) in 
this small open set must be part of this curve... but how small is small? Notice that 
the following two difficulties need to be taken care of: (i) The fact that 0(ti, ai, bi) 
is near ^ and 0 is not constant along the curve, does not imply that 0 eventu¬ 
ally reaches the value ^ on this curve, (ii) The fact that Q{ti,ai,bi) < ^ and 
0 (t 2 , 02 )^ 2 ) > ^ for two nearby points (ti,ai,6i) and (^ 2 , 025 ^ 2 ) that satisfy the 
equations R{t,a,b) = 0 = F{t,a,b), does not guarantee that 0 eventually reaches 
the value ^ due to the fact the point (ti,ai,6i) may no be in the same connected 


component of the curve of solutions that contains (t 2 , 02 ) ^ 2 )- Theorem 3.1 helps to 
solve these two difficulties. 


The idea of the proof of Theorem 6.1 is the following. We consider a box in the 
space, see Figure 6.1 We use Theorem 3.1 to show that every plane 6 = c in this 


box contains only one solution of the equation R{t,a,b) = 0 = F{t,a,b). We use 
the Poincare Miranda Theorem to show that there are reduced periodic solutions on 
three rectangles contained in planes of the form b = bo, b = bo — sb and b = bo + sb. 
We finish the prove by showing that 0 evaluated on the solution of the equations 
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Figure 6.1: Points {t,a,b) in the curve represent reduced periodic solutions of the 
three body problem. We proved the existence of a point in this curve in each one 
of the three small rectangles shown in the picture and we prove that at the point 
in the curve contained in the first rectangle the value of the function 0 is smaller 
than ^ and at the point in the curve contained in the last rectangle the value of the 
function 0 is bigger than ^. The proof of the existence of a point in the curve on the 
rectangle in the center was needed to implement our theorem regarding the Implicit 
Function Theorem, which, in this case, it guarantees the existence of a unique point 
on the curve in every plane b = c. 
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R{t, a, 6) = 0 = F(t, a, b) contained on the rectangle with b = bo — sb is less than ^ 
and 0 evaluated on the solution of the equations R{t, a, 6) = 0 = F{t, a, b) contained 
on the rectangle with b = bo + sb is greater than 


Now we are ready to state the four lemmas mentioned in the beginning of this section 
that lead to the proof of the main Theorem 6.1 The first lemma shows that there 


is a point in the curve on the box in Figure 6.1 in the rectangle in the center. 


Lemma 6.3. There is a solution of the equations R{t,a,b) = 0 = F{t,a,b) with 
b = bo, \a — uqI < da and \t — to\ < dt 


The second and third lemmas show that there is a point in the curve on the initial 
and final rectangles |6.1[ 

Lemma 6.4. For some {t, a, b) with b = bo — sb, |a — (ao + sa)| < da and \t — {to — 
st)| < dt, there is a solution of the equations R{t,a,b) = 0 = F{t,a,b) such that 
Q{t,a,b) < 

Lemma 6.5. For some {t, a, b) with b = bo + sb, |a — (ao — sa)| < da and \t — (to + 
st)\ < dt, there is a solution of the equations R{t,a,b) = 0 = F{t,a,b) such that 
e{t,a,b) > 

The fourth lemma guaranties the set for values of (t, a, b) given by 


[to — 6(st + dt),to + 6(st + dt)] X [ao — 3(sa + da),ao + 3(sa + da)] x [bo — sb, bo + sb] 

is small enough to only allow one connected curve on it as the solution for the 
equation R{t, a,b) = 0 = F{t, a, b). 

Lemma 6.6. For every b G [bo — sb, bo + sb] there exists a unique solution (t, a, b) 
of the equations R{t, a,b) = 0 = F{t, a, b) with |a — ao| < 3(sa + da) and |t — to| < 
6{st + dt). 

7 Bounds 

In this section we will define 3 differential equations that will help us find bounds 
for the functions F{t, a, b), R{t, a, b), Q{t, a, b) and their partial derivatives. We will 
be using the functions (fi defined in section 
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7.1 The vector field W 


Let us consider the vector field W = (x2, (/>i, X4, </)2, ^3)^, where the functions 
are those defined on section For any fixed a and b, the function 

t —>■ a, b), F{t, a, b), R{t, a, b),R{t, a, b), Q{t, a, 6)^ 

satisfies the differential equation Y{t) = W{Y(t)) with Y (0) = (0, 6,10, 0, 0)^ 

Proposition 7.1. Using the functions (pi defined in section\^ we have that the 
derivative matrix DW of the vector field W, and the derivative matrix DWi of the 
vector field Wi = DW W are given by, 


/ 

0 

1 

0 

0 

0 

\ 


/ 

(p4 

0 

4>5 

0 

0 

\ 


(p4 

0 

(pb 

0 

0 




(pw 

(p4 

4>12 

4>5 

0 



0 

0 

0 

1 

0 


and DWi = 


2(p^ 

0 

4>g 

0 

0 



2</>5 

0 

</>6 

0 

0 




2(/>12 

2<p5 

4>14 

4>6 

0 


V 

0 

0 

(p7 

0 

0 

/ 


V 

0 

0 

(pl5 

(p7 

0 

/ 


7.2 The vector field G 

Let us consider the vector held G = {x 2 , (pi, X 4 , (p 2 , xq, (piQ, xg, (pn, cpis)'^, where the 
function (pi are those dehned on section For any hxed a and b, the function that 
sends t to 


(^F{t, a, b),F{t, a, b), R{t, a, b), Rft, a, b),Fa(t, a, b),Fa{t, a, b),Ra{t, a, b), Ra{t, a, 6), 0a(t, a, b 
satishes the differential equation Y (t) = G{Y {t)) with Y (0) = (0, b, 10,0,0, 0, 0, 0, 0)^ 


Proposition 7.2. Using the functions (pi defined on section we have that the 
derivative matrix DG of the vector field G, and the derivative matrix DGi of the 
vector field Gi = DG G are given by, 
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and 


DGi 


DG 


/ 


V 



/ 

0 

1 

0 

0 

0 

0 

0 

0 




04 

0 

05 

0 

0 

0 

0 

0 

0 



0 

0 

0 

1 

0 

0 

0 

0 

0 



205 

0 

06 

0 

0 

0 

0 

0 

0 

= 


0 

0 

0 

0 

0 

1 

0 

0 

0 



019 

0 

020 

0 

04 

0 

05 

0 

0 



0 

0 

0 

0 

0 

0 

0 

1 

0 



2020 

0 

022 

0 

205 

0 

06 

0 

0 


V 

0 

0 

023 

0 

0 

0 

07 

0 

0 / 


04 

0 

05 

0 

0 

0 

0 

0 

0 

010 

04 

012 

05 

0 

0 

0 

0 

0 

205 

0 

06 

0 

0 

0 

0 

0 

0 

2012 

205 

014 

06 

0 

0 

0 

0 

0 

019 

0 

020 

0 

04 

0 

05 

0 

0 

030 

019 

034 

020 

010 

04 

012 

05 

0 

2020 

0 

022 

0 

205 

0 

06 

0 

0 

2034 

2020 

039 

022 

2012 

205 

014 

06 

0 

0 

0 

040 

023 

0 

0 

015 

07 

0 


\ 




7.3 The vector field U 

Let us consider the vector field 


U = {X2, 4'1,X4, (j)2, XlQ, (f>41,Xi2, 

where the function (/>* are those defined on section For any fixed a and b the 
function 


t —y 


(^F{t, a, b),F{t, a, b), R{t, a, b),R{t, a, b),Fb{t, a, b), Fb{t, a, b), 
Rb{t,a,b),Ri,{t,a,b)^ 


satisfies the differential equation Y{t) = U{Y{t)) with Y (0) = (0, b, 10, 0,0,1,0,0)^ 
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Proposition 7.3. Using the functions (fii defined in section we have that the 
derivative matrix DU of the vector field U, and the derivative matrix DUi of the 
vector field Ui = DU U are given by, 


and 


DUi 



( ° 

1 

0 

0 

0 

0 

0 

0 \ 


<i>A 

0 

(p5 

0 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 

0 

0 

DU = 

2(/)5 

0 

4>6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 


4>A3 

0 

4>44 

0 

04 

0 

05 

0 


0 

0 

0 

0 

0 

0 

0 

1 


V 2(?!)44 

0 

045 

0 

205 

0 

06 

0 / 



/ 04 

0 

05 

0 

0 

0 

0 



010 

04 

012 

05 

0 

0 

0 

0 


205 

0 

06 

0 

0 

0 

0 

0 


2012 

205 

014 

06 

0 

0 

0 

0 


043 

0 

044 

0 

04 

0 

05 

0 


048 

043 

050 

044 

010 

04 

012 

05 


2044 

0 

045 

0 

205 

0 

06 

0 


\ 2050 

2044 

053 

045 

2012 

205 

014 

06 / 


7.4 Reduced Periodic Solutions 


In this section we use theorem 2.3 and the Poincare -Miranda Theorem to prove the 


existence of three reduced periodic solutions, one on each rectangle in figure 6.1 


Remark 7.4. Several lemmas in this section will he using the Round Taylor method 


and therefore it will be using Theorem 2.3 to estimate the values of the solution of 
the ode’s. If we take a look a the hypothesis of this lemma we notice that there is a 
number e that has to be greater than M^h + H. This e will be in all proofs that 
use the Taylor method in this paper. Some of these lemmas contains the variable e, 
in each case, it just refers to a small number giving an estimate of the error. 

2677451 


For any a G [ao — 3(sa -|- da),ao + 3(sa + da)]. 


Lemma 7.5. Let e = tqooooooo- 
b G [bo — sb, Iq + sb] and t G [to — 6(st -|- dt),to + 6(st + dt)] we have that 


\Qa{t,a,b) - 


536760312951 

20000000000000 ^ 


< e 
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and, 


Fa{t,a,b) 


3032500537707 

10000000000000 


Fait, a, b) 


11824770099363 

25000000000000 


Ra{t,a,b) 


68073031375453 

25000000000000 


Rait,a,b) 


164497338366219 

100000000000000 


Proof. Let us consider the following intervals, 


1 

62' 

72 = 

1 

3' 

7^0 — 

'189 

1001 ' 

h = 

33 

1 

Jr: — 

1 

31 

100 

25 

100 

2 

^3 

20 

100 

100 ’ 

100 

-^5 

100 ’ 

100 


1 12' 

h — 

1 69' 

To — 

1 42' 

Jg — 

1 3 ■ 

~100’ ^ 

~ioo’ 

-^8 — 

~ioo’ 

~100’ 


A direct computation using the bounds in section shows that if 


42 305541 249309 _ 778131 _ 246743 

Moa — —, AToo — -, Kin — -, Flig —- ATjo —- 

25’ ^ 125000’ ^ 100000’ ^ 1000000 ® 250000 

374443 133409 1345793 2429239 

® 500000 ® 100000 ^ 1000000 ® 1000000 

833241 _ 3298559 _ 182893 __i4 , to 

^ 500000’ ^ 1000000’ 1000000’ ’ 30000 ’ 

then, for values of (xi,...,X8) with Xi G R we have that \GR < M^g where 
G^,... ,G^ are the entries of the vector field G; IG 2 I < Mjg where G 2 ,..., G® are 
the entries of the vector field G 2 = DGi G. Recall that Gi = DGG. Moreover we 
have that 


|77G| < K^g, \DGi\ < Kig 

The Round Taylor method of order 2 using the vector field G with a = ao, k = 30000 
and initial conditions T( 0 ) = ( 0 , 60 ) 10 , 0 , 0 , 0 , 0 , 0 , 0 )^ produces a sequence 
with Z30000 equal to 
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247458249564811 13245901 189061430242601 1795639 

100000000000000’ 100000000000000 ’ 20000000000000 ’ 12500000000000’ 
3032500537707 11824770099363 68073031375453 164497338366219 536760312951 

10000000000000’ 25000000000000 ’ 25000000000000’ 100000000000000’ 20000000000000 


In this case H = 


M 

6 + 


[^Lkh _ ~ 0.0000019 where Mg = yjM'^^ + • • • + M|^ 

and L = K^g + Kig^. A direct verification shows that, for every j = 1,... ,8 and 
., 30000, the entry of Zi is within a distance e of the boundary of the 


* = 1 , 

interval Ij. Also we have that the last entry of Zi is within a distance of the 
boundary of the interval Jg. By Theorem 2.3 we conclude that for alH = 1,..., /c, 


7/7 . 7/7 7/7 . oh oh . oh 

F{-, ao, 6 o), F(-, ao, bo), R{^,ao, bo), R{^,ao, bo), T,(-, ag, bo),Fa{-,ao, bo), 
K K K K rC K 

Ra{-j^,ao, bo), Rai-j^, do, bo), 0 a(^, Go, bo) 


is within a distance H of Notice that for any function p, under the assumption 
that \p{t, ao, 6 o)| < Mp for all r between t and to, we have that 


\p{t, a, b) - p{to, ao, bo)\ < \p{t, a, b) - p{t, ao, bo)\ + \p{t, ao, bo) - p{to, ao, 6 o)| 

< \p{t, a, b) - p{t, ao, bo)\ +Mp{t- to), 


We will use the observation above to finish the proof of the lemma. We will bound 
\p{t, a, b) — p{t, ao, 6 o)| using Theorem 2.4 with the vector fields G and Go, where Go 
is the vector field G with a replaced by ag. A direct computation shows that 


6 G — G — Go — ( 0 , 0 , 0 , (/> 54 , 0 , 0 , 0 , (^> 57 ) 

Using the information on section we obtain that Idd < gggggg- Therefore, using 
Theorem |2.4| we conclude that the values of the solution of the differential equation 
using G (with a general a and b) compare with those of the solution of the differential 
equation using Go differ by less than 


1^ _ 5 g|e^ 09 (io+ 6 (st+di)) _|_ /^Kog{to+6{st+dt)) 

500000Agg V 


267131 

10000000 
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Therefore we have that for any a G [ao—3(sa+(io), ao+3(sa+da)], h G [bo—sb, 6o+'S^] 
and t G [to — 6(st + dt),to + 6{st + dt)] 


Fa{t,a,b) 


3032500537707 , 267131 , , , _ ~ 2677451 

- <-h Gist + dt)Mn„ + H < - 

10000000000000' 10000000 ^ - jqoOOOOOO 


We have similar computations for the functions 0^, Fa, Ra and Ra- This finishes 
the proof. □ 

Corollary 7.6. For any a G [ao —3(sa + (ia), ao + 3(sa+da)], b G [bo — sb,bo + sb] and 
t £ [to — G{st + dt), to + 6(st + dt)] , Fa{t, a, b) > 0, Ra{t, a, b) > 0, and |0a(t, a, fo)| < 

27 

1000 

Lemma 7.7. Let e = igggggg^- For any a G [oq — 3(sa + da),ao + 3(so + da)], 
6 G [6o — sb, bo + sb] and t G [to — 6(st + dt), to + 6(st + dt)] we have that 


\Rb{t,a,b) 


88229751956717 

100000000000000 


|Ffe(t,a,6) 


50798112898451 

100000000000000 


Proof. Let us consider the following intervals. 


h 


1 62 

Too’ ^ 


d2 


1 3 

Too’ 2 


h 


189 1001 
100 


33 1 

Too’ Too 


1 89 

Too’ Too 

A direct computation using section [^ shows that if 


do = 


1 101 

Too’ 


dio = 


1 101 

2’ Too 


dll = 


1 81 

loo’ 100 


3 1226931 2557349 778131 _ 246743 

^2’ " 500000 ’ 1000000’ 1000000 250000 

374443 ,, 133409 ,, 765259 ,, 533571 

Mou — -, — -, Ads^ — -, Adgji —- 

500000’ 100000’ 500000’ 200000 

1000000 200000’ ’ 30000’ 

then, for values of (xi, X2, xg, X4, xg, xio, xn, X12) with x* G R we have that j < 
Mog where . ,U^ are the entries of the vector field U; jld^l < FIju where 
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U 2 , ■ ■ ■ ,U 2 are the entries of the vector field U 2 = DUi U. Recall that Ui = DU U. 
Moreover we have that 


\DU\ < Kou, \DUi\ < Kiu 

The Round Taylor method of order 2 using the vector field U with a = ao, k = 30000 
and initial conditions Y (0) = (0, bo, 10, 0, 0,1,0,0)^ produces a sequence with 

^30000 equal to 


247458249564811 13245901 189061430242601 1795639 

100000000000000 ’ 100000000000000 ’ 20000000000000 ’ 12500000000000’ 
25138479462137 50798112898451 20014508374143 88229751956717 \ 
12500000000000’ 100000000000000’ 25000000000000’ 100000000000000 J 


In this case H = 




/^Lkh — 1 "v ^ ^09 _ where A7 

Y ^ 100000000 wneie ivi-^ 


= + ^ 


+ Ml 


8u 


and L = Kqu + Kiu^- A direct verification shows that, for every j = 1,..., 8 and 
i = 1,, 30000, the entry of Zi is within a distance of the boundary of the 
interval Ij when j = 1, 2,3,4 and the entry of Zi is within a distance of the 
boundary of the interval Ij+A when j = 5,6, 7,8. By Theorem 2.3 we conclude that 
for all z = 1,..., /c. 


7/7 . 7/7 7/7 . oh oh . oh 

F{-,ao,bo),F{-,ao,bo),R{^,ao,bo),R{-,ao,bo),F,{-,ao,bo),Fh{-,ao,bo), 

K K K K K fv 

Rb{j^,ao, bo), Rb{^,ao, bo) 


is within a distance H of Zi. Notice that for any function p, under the assumption 
that \p{t, ao, 6o)| < Mp for all r between t and to, we have that 


\p{t, a, b) - p{to, ao, bo)\ < \p{t, a, b) - p{t, ao, bo)\ + \p{t, ao, bo) - p{to, ao, 6o)| 

< \p{t, a, b) - p{t, ao, 6o)| +Mp{t- to). 


We will use the observation above to finish the proof of the lemma. We will bound 
\p{t, a, b) — pit, ao, bo)\ using Theorem 2.4 with the vector fields U and Uq, where Uo 
is the vector field U with a replaced by oq. A direct computation shows that 


5U — U — Uo — (0,0,0, (/>54,0,0,0, (/iss) 
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Using the information on section we obtain that |(5[7| < 25000 • Therefore, using 
Theorem |2.4| we conclude that the values of the the solution of the differential equa¬ 
tion using U (with a general a and b) compare with those of the solution of the 
differential equation using Uq differ by less than 


|5 _ lj^^QKou(to+6(st+dt)) 


1 


25000ii:o„ 


^Kou{to+^{st+dt)) 


-1 < 


1281341 

50000000 


Therefore we have that for any a G [ao—3(sa+do), ao+3(sa-|-da)], b G [bo—sb, 6o+'S^] 
and t G [to — 6(st + dt), to + 6(st + dt)] 


\Fb{t,a,b) 


50798112898451, 1281341 , - 2568201 

- < -h 6(st + dt)Mnu -|- id < - 

100000000000000' 50000000 ^ - 100000000 


We have a similar computation for the function Rb- This finishes the proof. □ 

Lemma 7.8. Let e = 40000000 • Tor any a G [oo — 3(so -|- da),ao + 3(so -|- da)], 
6 G [60 — s 6 , 60 + st'] o.nd t G [to — 6 (st -|- dt), to + 6 (st + dt)] we have that 


\F{t,a,b) 


247458249564811 

100000000000000 


|d?(t, a, b) 


189061430242601 

20000000000000 


Proof. Let us consider the following intervals J 5 = [— 3 ^, 5 ] and di, I 2 , I 3 , h defined 


on the proof of Lemma 7.5 


A direct computation using the bounds on section [^ shows that if 


3 25282 260901 19453263 14977713 

^ 2 ’ ^ 15625’ 200000’ 25000000 20000000 

493485626283 _ 1334089805457 5396988231 

M2w — - = -, Mbw — - 1 

500000000000 1000000000000 125000000000’ 

then, for values of (xi,...,X 4 ) with Xi G R we have that jW-^l < Mq^ where 
W^, ..., are the entries of the vector field W ; | < Mj^, where W 2 , ■ ■ ■, W 2 are 

the entries of the vector field IT 2 = DWi W. Recall that lUi = DW W. Moreover 
we have that 


|dRU|<ddo^, |dRTi|<idi^ 


28 



















The Round Taylor method of order 2 using the vector field W with H = 10“^^, 
h = ^QQQQ , a = ao, k = 30000 and initial conditions Y (0) = (0, bo, 10, 0, 0)^ produces 
a sequence with Z 30000 equal to 


/247458249564811 


13245901 


189061430242601 


1795639 


12217304404331 


V100000000000000’ 100000000000000’ 20000000000000 ’ 12500000000000’ 10000000000000 


In this case H = 


6 ^ h 

L 
h 


(^^Lkh _ 1 ) < 


127 


1000000000 


where M,, 


= + ■■■ +Ml 


5w 


and L = Kow + Ri«) 5 • A direct verification shows that, for every j = 1,..., 4 and 

of the boundary of the 


i = 1,, 30000, the j*" entry of Zi is within a distance 
interval Ij. Also we have that the last entry of Zi is within a distance of the 


1000 


boundary of the interval J 5 . By Theorem 2.3 we conclude that for alH = 1,..., A: 


ih • ih ih ’ ih ?/? 

F{-,ao, bo),F{-, ao, bo), R{^,ao, bo),R{^,ao, bo), 0a(^, ao, bo) 

K rC K rv fv 


is within a distance H oi Zi. Notice that for any function p, under the assumption 
that \p{t, ao, 6 o)| < FIp for all r between t and to, we have that 


\p{t,a,b) - p(to,ao,^o)| < \p{t,a,b) - p{t,ao,bo)\ + \p{t,ao,bo) - p{to,ao,bo)\ 

< \p{t, a, b) - p{t, ao, 6 o)| +Mp{t- to). 


We will use the observation above to finish the proof of the lemma. We will bound 
\p{t,a,h) — p{t,ao,bo)\ using Theorem 2.4 with the vector fields W and Wo, where 
Wq is the vector field W with a replaced by oq. A direct computation shows that 


5W = W-Wo = (0,0,0,<^ 54 ^ 55 ) 

Using the information on section we obtain that \8w\ < iggggoo’ Therefore, us¬ 
ing Theorem |2.4| we conclude that the values of the the solution of the differential 
equation using W (with a general a and b) compare with those of the solution of the 
differential equation using Wo differ by less than 


\h — k \pKow(to+6ist+dt)) I _ yy _ ( Kow{to+S{^t+dt)) 

' ' lOOOOOORo^; V 


827737 

250000000 
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Therefore we have that for any a G [ao—3(sa+(io), ao+3(sa+da)], h G [bo—sb, 6o+'S^] 
and t G [to — 6 (st + dt),to + 6{st + dt)] 


\F{t,a,b) 


247458249564811, 827737 , -- 134567 

- <-h Q(st + dt)Mou] + H < - 

100000000000000' 250000000 ^ ; U«; _ 4 qqqqqqq 


We have a similar computation for the function R. This finishes the proof. □ 

Definition 2. We will denote by Zw{t,b,a,k) G the last vector in the sequence 
{^i}i=i produced using the Round Taylor Method of order 2 using the vector field 
W with H = 10“^“^, ^ ~ i initial conditions (0,6,10,0,0)^. We denote by 

Hw{t, k) = I'*' - l) where = y /H-h and L = Kqu, + 

Kiw^. If we need to use H = lO”*^ instead of H = 10“^^, then we will use the 
notation Zw{t, 6, a, k, q) G R^ and H\Y{t, k, q). 


Remark 7.9. For all the following Lemmas that use the Round Taylor Method, it 
can be directly verified that all the conditions of Theorem 2^ are satisfied; therefore 
Zy^{t,b,a,k) is within a distance H-\^{t,k) of 


{F{t, a, b),F{t, a, b), R{t, a, b),R{t, a, b), Q{t, a, 6))^. 

For the sake of keeping the flow of the proof, we will left out all the details in the 
proofs of these lemmas. 


Corollary 7.10. For any a ^ [oq — 3(sa + da), oq + 3(sa + da)], 6 g [6o — s6, 6o + s6] 


andtG [to — 6 (st + dt), to + 6 (st + dt)] we have that 


120689 

250000 


< 0 (t, a, 6 ) < 


483453 
1000000 ’ 


163169 

200000 


< F{t, a, b) < 


813693 

1000000 


and 


4593 

12500 


< R{t, a, b) < 


18653 

50000 


Proof. As in Lemma 


7.8 


let e = 


134567 


40000000 ' Notice that the values of F{t,a,b) are 
bounded by the minimum and maximum of the function — when 

2474582495648111 ^ 1890614302426011 


< e 


< e. 


■ (4x2+a;2)3/2 

100000000000000 1 ^ 1*3 — 2 lxwoool)booooo I ^ bound for F follows as a 

direct application of the Lagrange Multiplier method. Since the values of R{t, a, b) 
are bounded by the minimum and maximum of the function tllSS— ^ 


‘^3 ^3 (4a;2+x2)3/2 ’ 

we can use the same argument for R{t, a, b). Finally since the values of 0(t, a, 6 ) are 
bounded by the minimum and maximum of the function we can also use the 


same argument for 0 (t, a, 6 ). 


□ 
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Lemma 7.11. F{tQ — dt, a, bo) > 0 for all a G [oq — da, ao + da], F{to + dt, a, bo) < 0 
for all a £ [ao — da, oq + da], R{t, oq — da, bo) < 0 for all t G [to — dt, to + dt] and 
R{t, ao + da, bo) > 0 for all t £ [to — dt, to + dt]. 


Proof. By Corollary 7.6 we have that Fa{to—dt, a, bo) > 0 for all a £ [ao—da, oo+da]. 
Therefore, in order to show that F{to — dt, a, bo) > 0 for all a £ [oq — da, ao + da], it 
is enough to show that F{tQ — dt, ao — da, bo) > 0. A direct computation shows that 
^w(io — dt, ao — da, bo, 35000) equals to 


7733069351623 25787091 189061375789453 44841643 30543236182739 \ 

3125000000000’ 20000000000000’ 20000000000000 ’ ~ 20000000000000 ’ 25000000000000 J 


and 


Therefore, 


Fw{to — dt, 35000) < 


94851 

1000000000000 


F{to — dt, ao — da, bo) > 


25787091 


94851 


20000000000000 1000000000000 


> 0 


As pointed out above, we conclude that F{to — dt, a, bo) > 0 for all a G [ao — da, ao + 
da]. In the same way we have that F{to + dt, a, bo) < 0 for all a G [ao — da, oq + da] 
because a direct computation shows that Z\Y{to + dt, ao + da, bo, 120000) equals to 


1546614123963 198811 945307243792047 16995193 61086532113189\ 

625000000000 ’ ~6250000000000’100000000000000’20000000000000’500000000000007 


and 


Hw{to + dt, 120000) < 


5651 

200000000000 


Therefore, 


, , , N 198811 

Fito + dt, ao + da, bo) < - 

u -r , u w , u; 625000000000C 
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Let us show that R{t, a^ — da, bo) < 0 for all t G [tQ—dt, to+dt], A direct computation 
shows that Zw{to + dt, oq — da, bo, 35000) equals to 


/123729109625969 196972647 945306878946787 19027313 122173137972697 

V 50000000000000 ’”100000000000000’100000000000000’”25000000000000’100000000000000 


and 


Hw{to + dt, 35000) < 


94851 


1000000000000 


Therefore, 


R{to + dt, 


0-0 


, , , 19027313 94851 

- da, bo) < -1-< 0 

^ 25000000000000 1000000000000 


Since, R{t, ao — da, bo) > 0 (see Corollary 7.10), we conclude that R{t, ao — da, bo) < 0 
for all t G [to ~ dt,to + dt]. Let us show that R{t,ao + da, bo) > 0 for all t G 
[to — dt, to + dt]. A direct computation shows that Zy/{to — dt, ao + da, 6o) 35000) 
equals to 


24745827990179 55883369 189061484706171 52393253 61086475049353\ 

10000000000000’ 25000000000000 ’ 20000000000000 ’ 50000000000000’ 50000000000000 ) 


and 


Hw{to ~ dt, 35000) < 


94851 

1000000000000 


Therefore, 


i?(to — dt, ao + da, bo) < 


52393253 


94851 


50000000000000 1000000000000 


> 0 


Since, R{t, ao+da, bo) > 0 (see Corollary 7.10), we conclude that R{t, ao+da, bo) > 0 
for all t G [to — dt, to + dt]. This hnishes het proof of the Lemma. 


□ 


Remark 7.12. As a Corollary, using the Poincare-Miranda theorem, we obtain a 
proof of Lemma |d.3t 
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Lemma 7.13. F{to — st — dt, a, bo — sb) > 0 for all a £ [oq + sa — da, oq + so + da], 
F{to — st + dt, a, bo — sb) < 0 for all a £ [oq + so — da, oq + so + da], R{t, oq + so — 
da, bo — sb) < 0 for all t £ [to — st — dt, to — st + dt] and R{t, oq + so + da, bo — sb) > 0 
for all t £ [to — st — dt, to — st + dt]. Moreover, for all all t £ [to — st — dt, to — st + dt] 
and all a £ [ao + sa — da, oq + so + da], 0(t, a, bo — sb) < ^. 


Proof. The proof is similar to that of Lemma 7.11 In this case F{to — st — dt, a, bo — 
sb) > 0 for all a £ [oq + so — da, oo + so + do] because Z]v{to — st — dt, oo + so — 
do, bo — sb, 35000) equals to 


247454580109467 14064311 945308716843341 27791851 7635805749249 \ 

100000000000000’ 50000000000000’ 100000000000000’ ”50000000000000’ 6250000000000 ) 


and 


Hw{to ~ st 


dt, 35000) < 


94849 

1000000000000 


Therefore, 




F{to — st + dt, a, bo — sb) <0 for all a £ [oq + sa — da, oq + so + do], because 
Zw{to — st + dt, oo + sa + da, bo — sb, 120000 ) equals to 


49490920135587 4841111 29540903186787 21447973 122172932415207\ 

20000000000000’~100000000000000’ 3125000000000 ’ 25000000000000’100000000000000 ) 


and 


Hw{to 


st + dt, 120000 ) < 


5651 

200000000000 


Therefore, 


F{to — st + dt, ao + sa + da, 




hr. . 
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R{t, ao + sa + da, bo — sb) > 0 for all t £ [to — st — dt, to — st + dt], because Zw{to — 
st — dt, ao + so + da, bo — sb, 35000) equals to 


123727300365019 12058321 189061780400321 11254969 61086446906439\ 

50000000000000 ’ 20000000000000 ’ 20000000000000 ’ 20000000000000 ’ 50000000000000 ) 


and 


Therefore, 


Hw{to — st — dt, 35000) < 


94849 

1000000000000 


R{to — st — dt, ao + sa + da, bo — sb) > 


11254969 


94849 


R{t, ao + sa — da, bo — sb) < 0 for all t £ [to — st — dt, to — st + dt], because Zw{to — 
st + dt, ao + sa — da, bo — sb, 35000) equals to 


4949091602187 18526257 945308716843239 25964629 122172930636361 \ 

2000000000000’”50000000000000’ 100000000000000’”100000000000000’ 100000000000000 ) 


and 


Therefore, 


Hw{to — st — dt, 35000) < 


1897 

20000000000 


R(to — 


In order to show that &{t,a,bo — sb) < we first noticed that Zw{to — st,ao + 
sa, bo — sb, 35000) equals to 


247454590419723 1161959 945308809422539 473597 122172912224601 \ 

100000000000000 ’ 10000000000000 ’ 100000000000000 ’ 3125000000000’ 100000000000000 ) 
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and 


Hw{to ~ st, 35000) < 


94849 

1000000000000 


Therefore, 


77 r 122172912224601 94849 Tvr 

0{to — st,an + sa,br) — so) -<-1- 

^ ’ ’ ^18 100000000000000 1000000000000 18 

122 1 7292 1 709501 77r 
100000000000000 ~ 18 


Using the mean value theorem, and the bounds |0a| < 


483453 

1000000 


and | 0 | < found 


in Corollaries 7.10 and 7.6, we obtain that for any t G [to — st — dt, to — st + dt] and 
a G [ao + sa — da, ao + sa + da] 


0 (t, a, bo — sb) — — = 0 (t, a, bo — sb) — 0 (t, ao — sa, bo — sb) + 

18 

0 (t, ao + sa, bo — sb) — 0 (to, ao + sa, bo — sb) + 

Ytt 

0 (to, ao + sa, bo - sb) - — 

lo 

483453 27 122172921709501 77r 

1000000 ““ 1000 100000000000000 ~ 18 

< 0 


□ 


Remark 7.14. As a Corollary, using the Poincare-Miranda theorem, we obtain a 
proof of Lemma Id.f] 

Lemma 7.15. F(to + st — dt, a, bo + sb) > 0 for all a G [ao — sa — da, ao — sa + da], 
F{to + st + dt, a, bo + sb) < 0 for all a G [ao — sa — da, ao — sa + da], R{t, ao — sa — 
da, bo + s 6 ) < 0 for all t G [to + st — dt, to + st + dt] and R{t, ao — so + da, bo -hsb) > 0 
for all t G [to + st — dt, to + st + dt]. Moreover, for all all t G [to + st — dt, to + st + dt] 
and all a G [ao — sa — da, ao — sa + da], 0(t, a, 6 o + s 6 ) > ^. 


Proof. The proof is similar to that of Lemma 7.13, In this case F{to + st - 
sb) > 0 for all a G [ao — so — da, ao — so + da] because Zw{to + st — dt 
da, bo + sb, 35000) equals to 


- dt, a, bo + 
, oo — so — 
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247461898442221 1573331 189061080098229 57188327 15271644451087 

100000000000000’5000000000000’ 20000000000000 ’”100000000000000’12500000000000 


and 


Hw{to + st 


dt, 35000) < 


23713 

250000000000 


Therefore, 


F{to + st + dt,a,bo + 6) < 0 for all a € [oq + so — da,ao + so + da], because 
Zw{to + st + dt, oo — so + da, bo + sb, 120000,15) equals to 


( 


1237309595658901 904301 9453055857540481 


500000000000000 ’ 62500000000000’ 1000000000000000’ 

842419313 244346392152949 , 


1000000000000000 ’ 200000000000000 


and 


Therefore, 


Hwito + st + dt, 120000,15) < 


77 

8000000000 


R{t, ao — sa + da, bo + sb) > 0 for all t £ [to + st — dt, to + st + dt], because Zw{to + 
st — dt, ao — sa + da, bo + sb, 35000) equals to 


49492383812687 63630339 945305585649929 1366741 24434631486741 \ 

20000000000000’ 100000000000000’100000000000000’2500000000000’20000000000000 ) 
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and 


Therefore, 


Hw{to + st 


dt, 35000) < 


23713 

250000000000 


R{to — st — dt, oo + so + da, bo — sb) > 


1366741 


23713 


2500000000000 250000000000 


> 0 


R{t, ao + sa — da, bo — sb) < 0 for all t £ [to + si — dt, to + st + dt], because Zw{to + 
st + dt, ao — sa — da, bo + sb, 35000) equals to 


123730949221087 33715619 945305400490959 6891931 1908956160269 \ 

50000000000000 ’”100000000000000’ 100000000000000’”25000000000000’ 1562500000000 y 


and 


Hw{to + st + dt, 35000) < 


23713 

250000000000 


Therefore, 


R(tn 


In order to show that Q{t,a,bo + sb) > we first noticed that Zw{to + st,ao — 
sa, bo + sb, 35000) equals to 


123730954376411 3739331 945305493070449 1355097 24434635169083 \ 

50000000000000 ’ 25000000000000 ’ 100000000000000 ’ 10000000000000 ’ 20000000000000 ) 


and 


Therefore, 


Hw{to + st, 35000) < 


23713 

250000000000’ 
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77r 24434635169083 23713 

0(tn — an + sa, On — so)->- 

^ ’ ’ ^18 20000000000000 250000000000 

24434637066123 lix 
20000000000000 ~ 18 


I'K 

18 


Using the mean value theorem, and the bounds |0a| < i^o^oo^q^o^q 
on corollaries 7.10 and 7.6, we obtain that for any t G [tg + st 
a G [oq — sa — da, og — sa + da] 


and | 0 | < found 
— dt, to + st + dt] and 


0 (t, a, 6 g — sb) 


— = &{t,a,bo - sb) - Q{t,ao - sa,bo - sb) + 

io 

0(t, ag + sa, bo — sb) — 0(tg, ag + sa, bo — sb) + 

hr 

0(tg, ag + SO, bo - sb) - — 

lo 

24434637066123 hr 483453 27 

^ 20000000000000 ~ 18 ~ 1000000 ““ ~ looo ® 
> 0 


□ 


Lemma 7.16. Let Si = 


8159 ^ _ 2677 r _ 16181 ? _ 8359 , _ 4549 
10000’ ~ 5000’ ^2 — igggg; ^2 — 5 ggg, t2 — 5 ggg, 


ei = 6 (st + dt), €2 = 3(sa + da), /xi = dt and ^2 = da. If pi 


nr , = ^ 2-/^2 

mi ^ m 2 ’ 


where mi = slsfZHl "^2 = ^ ^ [^o - eiUo + h], 

a G [og — 62 , ag + 62 ] and b G [bo — sb, bO + s 6 ] we have that 


61 <\F{at,a,b)\ < 5 i, \Fa{t,a,b)\ < ei, \Fb{t,a,b)\ < ei 


and 


\R{t,a,b)\ < € 2 , 82 < \Ra{t,a,b)\ < 82 , \Rb{t, a, b)\ < 62 

Moreover, sb < pi and sb < p 2 . 


Proof. This is a direct computation. It uses Lemmas 7.5 7.7 and 7.10 


□ 


Remark 7.17. The previous lemma, along with Theorem 3.1 implies Lemma 6.6 
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Remark 7.18. Since we have shown the existence of a curve of initial conditions 
that provides reduced periodic functions, then we have not only shown the existence 
of the periodic solution with 0{to) = ^ but we have shown the existence of infinitely 
many periodic solutions due to the fact that the function 6 is not constant along this 
curve, and, on any open interval, there are infinitely many numbers of the form ^ 
with p and q whole numbers. 

8 Dictionary of functions 

In this section we define the functions involved in the definition of the differential 
equations that we considered in this paper. The domain of the variable xi will be 
an interval containing possible values of the function F(t,a,b). In the same way, 
the variable X 2 is related with the function F{t, a, b); the variable X 3 is related with 
the function R{t,a,b); the variable X4 is related with the function R{t,a,b); the 
variable X 5 is related with the function Fa{t,a,b); the variable x% is related with 
the function Fa{t,a,b); the variable xj is related with the function Ra{t,a,b); the 
variable xg is related with the function Ra{t,a,b); the variable xg is related with 
the function Fh{t,a,b)-, the variable xio is related with the function Fh{t,a,b)] the 
variable xn is related with the function Ri,{t,a,b); the variable X 12 is related with 
the function Ri,{t,a,b). As in the previous section, s = \J4xf + Xg. For each one 
of these function {lbi,ubi} where Ib^ is a lower bound and ub* is a lower 

bound for the function fii = (/)j(xi, X2, X3, X4, X5, xg, X7, xg, xg, xio, xn, X12, a) on the 
set V, where 

V = h X I2 X I3 X I4 X I5 X Iq X R X Ig X Ig X ho X hi x h2 x /13 

with 


j _\ 1 621 ^ _ r 1 3 j ^ _ ri 89 looij ^ _ r 33 1 1 ^ [ 1 31' 

^“[~ioo’2j ^ 100 J “ [~ioo’looj ^“["iTO’lOO 


r 1 12] r 1 69] ^ r 1 42] r 1 101] ^ ri lor 

100 25J [ 100 25J [ 100 25J [ 100 50 J [2’100 

_r 1 81 1 _ r 1 89 1 _ r43170106052787 43170844652787' 

““["lOO’iooJ “ [~100’lOoJ “ [ 10000000000000’ 10000000000000 
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We would like to point out that /13 = [oq — 3(sa + da),ao + 3(sa + da)]. 

Remark 8.1. In order to obtain each bound B{(pi) = {Ibi, ubi}, we first compute the 
minimum and maximum of the function, mi and ma respectively. Then we define 
Ibi = ^[W^mi] and ubi = , where k is a positive integer. Usually k is 

between 4 o,nd 12 and it is chosen according to the precision that we want for the 
bound. The reason for using Ibi, ubi and not the minimum and maximum is that ma 
and mi may have complicate expression in term of radicals and roots of polynomials. 

Remark 8.2. The bounds for each of the functions for which we are writing an 
expression in terms of the variables xi, ..., X 12 have been obtained by directly using 
the Lagrange multipliers method. For those functions (pi that we are written in 
term of the previous (pj and some Xk, we used the method of Lagrange multiplier by 
changing each (pj by a variable Uj and we used the fact that Uj has bound given by 
B{(Pj). 


8.1 


400x1 


n — 


Entry 2 of the vector field W. B(<pi)= {- } 25 ooo > SUO^} 


8.2 



25 200x3 

x| 


Entry 4 of the vector field W. B(<p 2 )= {" iSoV 1000000 } 


8.3 



Entry 5 of the vector field W. B(fi 3 )= 1 qooooo } 


8.4 



Entry (2,1) of the matrix dW. B{(p 4 )= - Soo } 
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8.5 


1200a;ia;3 


Entry ( 2 , 3 ) of the matrix dW. B{(j)5)= {“2^^, 50^} 


8.6 


= 50 


6 a^ 


Xn 


+ 


12 x1 




Entry ( 4 , 3 ) of the matrix dW. B((^6)= {“ 


83881 


108213 


200000 ’ 1000000 . 


8.7 07 = 


20a 


Xn 


Entry ( 5 , 3 ) of the matrix dW. 8(^7)= {-3^^,- 1000000 } 


8.8 


1200 {IQxlx^ - xi) 


Part of the entry ( 2 , 1 ) of the matrix dWi. B( 08 )= { — 500000 ’ 12^00^0 ) 


, 4800 (8a;? — 3a;ia;o) 

8.9 09 =-^ 


Part of the entry ( 2 , 1 ) of the matrix dWi. B( 09 )= {- 1000000 ’ 1000000 ) 


8.10 010 = a ;408 + a ;209 

Entry (2,1) of the matrix dWi. B( 0 io)= {- 1000000 ’ ) 


8.11 


011 = 


4800 (a;f — a;ia;|) 


Part of the entry ( 2 , 3 ) of the matrix dWi. B( 0 ii)= losi^} 
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8.12 012 = a;40ii + X2(j)8 


Entry ( 2 , 3 ) of the matrix dWi. B(0i2)= {- iqqoooo ^ } 


8.13 013 


12000^ _ 3000x1 1800x3 _ 150 

q'J / y »4 

o o ^3 


Part of the Entry ( 4 , 3 ) of the matrix dWi. B(0i3)= 


8.14 014 — X 4013 + 2 X 2011 


Entry ( 4 , 3 ) of the matrix dWi. B(0i4)= {- 1ooooom soMto} 


8.15 015 


60ax 


Xq 


4 


Entry ( 5 , 3 ) of the matrix dWi. B(0i5)= 4^} 


8.16 016 


400 (8x5X1 + 3X3X7X1 — X3X5) 


Entry 6 of the vector field G. B( 0 i 6 )= 1000000 } 


8.17 017 — —1007 4 “ 2 X 505 4 “ 2^706 


Entry 8 of the vector field G. B(0i7)= {-50^, ^qooooo } 


8.18 


10 (x3 — 2 axj ) 


X 


3 

3 


Entry 9 of the vector field G. B( 0 i 8 )= 1000000 } 


8.19 019 — X 509 + X 708 


Entry (6, 1 ) of the matrix dG. B{(j)ig)= {-5^^, Soo } 


9611 1 
40000 / 
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8.20 020 — 3^508 + 2^7011 


Entry ( 6 , 3 ) of the matrix dG. B(02o)= {- looooom 20U^} 


8.21 


21 — 


600a 




Part of the entry ( 8 , 3 ) of the matrix dG. B(02i)= - 1 oooqqo } 


8.22 022 — 021 + 2X501 + X7013 


Entry ( 8 , 3 ) of the matrix dG. B(022)= Hgg} 


r 961 32973 1 

\ 40000’ 500000 / 


8.24 024 = X808 + X6(j)g 


8.23 023 — 


60ax 


7 


Xn 


20 


^3 ^3 

Entry ( 9 , 3 ) of the matrix dG. B( 023 )= 


Part of the entry (6, 1 ) of the matrix dGi. B( 024 )= {~ 125000 ’ 2^ooroo } 


8.25 4800(16xf-27xix; + xS 

Part of the entry (6, 1 ) of the matrix dGi. B( 025 )= {- 25000m 
, 24000(16x?X 3 — 3xiXo) 

8.26 026 = -^-V-^ 

Part of the entry (6, 1 ) of the matrix dGi. B( 026 )= {- 1000000 ’ 1000000 } 

8.27 027 = X4025 + X2026 

Part of the entry (6, 1 ) of the matrix dGi. B( 027 )= { — 2^5{)roo ’ 200000 } 
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4800 (128x1 - 96x?x? + 3x1) 
8.28 028 = -^^ 


Part of the entry ( 6 ,1) of the matrix dGi. B(028)= {- loooooo ’ 5 ^} 


8.29 029 — 3^2 028 + X4026 


Part of the entry ( 6 ,1) of the matrix dGi. B(029)= {- 12 ^ 5000 , 1000000 } 


8.30 030 — 024 + a:7027 + a:5029 


Entry ( 6 ,1) of the matrix dGi. B(03o)= 1 ooS } 


8.31 031 — X608 + 3^8011 


Part of the entry (6,3) of the matrix dGi. B(03i)= {- looorom 1000000 } 


8.32 032 


24000 (Sx^xs — X 1 X 3 ) 


Part of the entry (6,3) of the matrix dGi. B(032)= {- 1000000 ’ 2 U^}- 


8.33 033 — X2025 + X4032 


Part of the entry (6,3) of the matrix dGi. B(033)= {- ioqoooo ’ 1000000 } 


8.34 034 — 031 + X7033 + X5027 


Entry (6,3) of the matrix dGi. B(034)= {- 5 ^, 


8.35 035 


2400ax4 


Part of the entry (8,3) of the matrix dGi. B(035)= {- 1000000 ’ sMo) 
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8.36 036 — 035 + 2^8013 + 2X6011 

Part of the entry (8,3) of dGi. B(036)= {-^^ 5 , 


8.37 


037 


600 


/ lOa^ 8x3 


108X3X3 12X3 1 \ 

n .7 I 


Part of the entry (8,3) of the matrix dGi. B(037)= { — 


8.38 038 — X4037 + 2X2032 

Part of the entry (8,3) of the matrix dGi. B(038)= {- 3000000 ^ loooroo } 


8.39 039 — 036 + 2X5033 + 2^7038 

Entry (8,3) of the matrix dGi. B(039)= {- 300000 V 1000000 } 


60 (x 3 (axg + X 4 ) — 40 X 4 X 7 ) 

040 = - 5 - 

x§ 

(9,3) of the matrix dGi. B( 04 o)= {- loVoro ^ 

400 ( 8 x 9 X 3 + 8 x 3 X 33 X 3 — x^g) 

<^41= 

6 of the vector field G. B( 043 )= 

8.42 042 = 2X905 + 2^1106 

Entry 8 of the vector field G. B(042)= lu^} 

8.43 043 = X1108 + X909 

Entry (6,1) of the matrix dG. B(043)= {- 30 ^) 000 ’ 


8.40 

Entry 

8.41 

Entry 
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8.44 044 = a;ii0ii + xg(j)s 

Entry (6,3) of the matrix dU. 500000} 

8.45 045 = 2X9011 + X11013 

Entry (8,3) of the matrix dU. B(045 )= ioooqqo } 

8.46 046 = Xi208 + a:io09 

Part of the entry (6,1) of the vector field dUi. B(046 )= {~ iootooq 1 31250 } 

8.47 047 = X4026 + a;2028 

Part of the entry (6,1) of the matrix dUi. B(047)= 1000000 ) 

8.48 048 = X11027 + X9047 + 046 

Entry (6,1) of the matrix dUi. B(048)= {-f^oM’ HM} 

8.49 049 = Xi20ii + a;io08 

Part of the entry (6,3) of the matrix dUi. B(049)= {- lopoooo ’ 200000 ) 


8.50 050 — Xg(p27 + a;ii033 + 049 

Entry ( 6 , 3 ) of the matrix dUi. B( 05 o)= 1000000 ) 


12000 ^ 


1200 x| 


2400x^^3 288 OOX 4 X 3 


Part of the entry ( 8 , 3 ) of the matrix dUi. B( 05 i)= { 


500000 ’ 40000. 


8.52 052 — 2xio0ii + a;i205i 

Part of the entry ( 8 , 3 ) of the matrix dUi. B( 052 )= {- iroooom 2^)000 ) 
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8.53 053 — 2X9033 + 2^11038 + 052 

Entry (8,3) of the matrix dUi. B(053)= {- Soq , Soq } 


8.54 054 = 


100(fl — Oo)(o “1“ Oq) 


Xn 


Entry 4 of VE - VEo- B(054)= {- 


19 19 


500000 ’ 500000. 


8.55 055 = 


10(a — Oq) 




Entry 5 of fE - VEo- B(055)= {- 50 ^, 20 ^} 


8.56 056 = 


100(a — Oq) ( 3 ax 7 + 809X7 — 2x3) 

xi 


Entry 8 of G - Gq. B(056)= 40 ^} 


8.57 057 = 


20 X 7(0 — Oq) 


Xq 


Entry 9 of G — Gq. B(057)= { — 


1000000 ’ 1000000 . 


8.58 


300 Xii(o — Oo)(o + Oq) 


Xq 


Entry 8 of [/ Uq- B(057)— { 1000000 ’ 1000000 } 
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